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Abstract a composite radar scattering model composed of the atmosphere radiative transfer model, 
and the ocean surface Bragg wave theory is developed to analyze the impact of hurricane rain on the nor- 
malized radar-backscatter cross section (NRCS) measured in the VV and cross-polarized C-band Synthetic 
Aperture Radar (SAR) channels. The model results are validated against SAR and SFMR measured wind 
speeds and rain rates for two hurricane cases. The contribution of rain to the NRCS is backscatter from two 
parts: the atmosphere column and the ocean surface. In the atmosphere, microwave attenuation and the 
rain-induced volume backscattering are simulated by the model. We find that the impact of raindrops in 
the atmosphere is almost negligible for the VV polarization, but important for the cross polarization. On the 
ocean surface, comparisons between our model and other existing models without rain lead to the conclu- 
sion that the VV polarization NRCS can be simulated reasonably well without considering the non-Bragg 
scattering mechanisms. Similar to the wave breaking mechanism, the microwave diffraction on the craters, 
crowns, and stalks, produced by rain drops, is also negligible for VV polarization. However, the non-Bragg 
scattering is important for the cross-polarized NRCS simulations. Finally, we performed simulations to under- 
stand the VV-polarized NRCS behavior under different wind speeds at various rain rates. 


1. Introduction 


Hurricanes or typhoons (all storms are referred to as hurricanes hereafter) have been frequently observed in 
spaceborne Synthetic Aperture Radar (SAR) images since the first SAR image became available in 1978 [Fu 
and Holt, 1982]. With the advantage of its high spatial resolution, relative large spatial coverage (especially 
with ScanSAR technology) and the capability of imaging the imprint of hurricanes on the ocean surface 
under almost all-weather conditions, the SAR instrument is suitable for detailed hurricane observations over 
the sea surface. Therefore, SAR images were adopted by many scientists to understand hurricane eye char- 
acteristics [Du and Vachon, 2003; Jin et al., 2014], morphology [Friedman and Li, 2000; Li et al., 2013], tracks 
[Zheng et al., 2015], intensity [Reppucci et al., 2010; Zhang et al., 2014], winds [Zhang and Perrie, 2012; Zhou 
et al., 2013; Li, 2015], and swell waves [Li et al., 2002]. However, rainfall associated with hurricanes was not 
addressed quantitatively in these studies. In general, rain has two significant effects on the normalized 
radar-backscatter cross section (NRCS) of the SAR: (1) raindrops induce volumetric scattering and attenua- 
tion in the atmosphere and (2) rain alters the roughness of the ocean surface. The aim of this study is to 
understand each of these mechanisms under hurricane conditions. 


On the ocean surface, a NRCS model including the rain damping effect on the wind-induced waves, and the 
rain enhanced ring wave spectra, was developed by Contreras and Plant [2006]. They added the two rain 
effects to a semiempirical wind wave spectrum model developed by Kudryavtsev et al. [2003] (hereafter 
denoted the KHCC wave spectral model). However, the sea surface waves that satisfy Bragg resonance con- 
dition are tilted by the long waves, so the local incident angles and polarizations are altered [Valenzuela, 
1968, 1978; Bass et al., 1968; Plant, 1990]. This tilting effect is merely dealt with by averaging over the scales 
of the long waves in the Bragg scattering model adopted by Kudryavtsev et al. [2003] and Contreras and 
Plant [2006], instead of using a probability density function (PDF) methodology. As the geometric coeffi- 
cients are zeros for the cross polarizations, the averaging process over the scales of long waves cannot 
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simulate the cross-polarized NRCS caused by Bragg wave resonance. Recently, Kudryavtsev et al. [2014] 
simulated the NRCS due to the Bragg resonance by adopting the PDF methodology. In the atmosphere, the 
raindrops attenuate the radar signal, and also intensify the NRCS through volumetric scattering. These two 
mechanisms were simply modeled by an empirical approach by Nie and Long [2007] without rain effects on 
the ocean surface. Thereafter, a more sophisticated physics-based radiative transfer model was developed 
by Xu et al. [2015], which is shown to give comparable results to those of the simple model of Nie and Long 
[2007]. Although Xu et al. [2015] simply modeled ring waves generated by rain with some success there is 
need for improvements; the mechanisms by which rain affects the ocean surface need further study. In Xu 
et al. [2015], the rain-induced NRCS damping was not included and the ring waves’ contribution was linearly 
added as a modulation to the empirical geophysical model function (GMF) without considering the actual 
physical progresses. In this study, rain’s effects on the ocean surface including the rain-induced ring waves 
and the damping due to the wind waves were modeled by the Bragg scattering theory and the semiempiri- 
cal wave spectrum. Thus, we can analyze each NRCS modulation mechanism quantitatively. 


In this paper, we investigate the hurricane rain effect on C-band SAR observations in terms of mechanisms 
related to the ocean surface, as well as those related to the atmosphere. Our approach is to formulate a 
composite model to include the rain effects in the atmosphere and also, the rain effects on the ocean sur- 
face. Section 2 describes our composite model and the rationale regarding why each piece is chosen. In sec- 
tion 3, the model is validated using the observed rain rates, wind speeds, and C-band VV-polarization 
NRCSs over two hurricanes. Finally, analysis of each mechanism effect on the SAR hurricane observation 
and conclusions are given in sections 4 and 5, respectively. 


2. The Radar Scattering Model 


In this study, we build a semiempirical radar scattering model that takes into account the impacts of both 
sea surface wind and rain on the NRCS. The model considers the total NRCS as the summation of backscat- 
ters from both the ocean surface and the atmosphere. As shown in Figure 1, there are five modules that are 
needed to account for these processes: (1) the composite sea surface Bragg model given by Valenzuela 
[1978] and Plant [1990] for its ability to simulate quad-polarized data; (2) the wind-driven KHCC gravity wave 
spectrum developed by Kudryavtsev et al. [2003] without the rain effect; (3) the rain damping effect on the 
KHCC wave spectrum given by Contreras and Plant [2006] and Tsimplis [1992]; (4) the ring wave spectrum 
developed by Le Méhauté [1988]; (5) an additive model developed by Nie and Long [2007] for the atmos- 
pheric part. With this new composite model, the mechanisms of the rain effect on the NRCS are analyzed. 


The composite Bragg model we formulated explains C-band NRCS measurements better and is less sensi- 
tive to the choice of the roughness spectral model than is the case for Ku-band measurements [Hwang and 
Plant, 2010]. Although other wind-wave spectra are probably more commonly used [Elfouhaily et al., 1997], 
the KHCC spectrum has an advantage. The main reason for choosing KHCC is that it is developed from a 
balance of wind input, parasitic wave generation, and wave dissipation, and therefore is convenient to mod- 
ify to include the rain damping effect on the wave dissipation. 


2.1. The Surface Scattering Model 

The Bragg resonance theory has been developed since the 1960s [e.g., Wright, 1966; Valenzuela, 1968; Bass 
et al., 1968; Valenzuela, 1978; Plant, 1990; Hwang et al., 2010]. The backscatter cross-section solutions for a 
slightly rough patch tilted by the longer waves are as follows: 
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where k is the radar wave number, y and ô are the tilting angles of the long wave in the radar incident plane 
and perpendicular to this plane, respectively, and 6; is the local incidence angle of the slightly rough patch: 
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Figure 1. Outline of the model. The red abbreviations mean the references: KHCC for Kudryavtsev et al. [2003], CP06 for Contreras and Plant 
[2006], T92 for Tsimplis [1992], L88 for Le Méhauté [1988], V78 for Valenzuela [1978], P90 for Plant [1990], and NLO7 for Nie and Long [2007]. 
6;=arccos [cos (0+1)cos ô] 
aj=sin 0; (2) 
x=sin (0+4), y=cos (0+4) 
where «, «, and y are parameters to simplify the formula, W(Kg,, Kgy) is the 2D-wave number variance spec- 
trum, gyy and g;, are the Bragg scattering geometric coefficients for VV and HH polarizations: 


(er—1) [e-(1+sin? 0) —sin? 6] 


Iw (4)= 
í [ecos 0+ e — sin? o) 
Imu (8)= (eV) (3) 


[cos 0+ vV &— sin? J f 


and where ¢, is the relative dielectric constant of seawater. Accounting for all the surface tilting, the back- 
scatter cross section per unit area of the sea surface is 


co(0) -| i co(0i)p(tan w, tan ô)d(tan y)d(tan ô) (4) 


where p(tany, tanò) is the joint probability density function (PDF) of the ocean surface slopes. We adopt the 
Gram-Charlier distribution [Cox and Munk, 1954] for the PDF function for p: 
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where 7 =tanw/s, and ¢ = tanô/se are the normalized upwind and crosswind slope components, 
respectively, and s? and sł are the crosswind and upwind Mean Square Slope (MSS) variables [Hwang 
et al., 2010]. The skewness also given by Hwang et al. [2010], as contained in the coefficients C21, C4o, 
C22, Coa, and Coz, increases with wind speed from nearly zero at low wind to c24 = —0.11, cy = 04, 
C22 = 0.1, Cog = 0.2, and Co3 = —0.42 at 14 m/s; a linear approximation is used in the implementation 
such that: 


C21) =—0.11U 49/14 (Uio < 14m/s) 
Co3 = —0.42U 19/14 (Uio < 14m/s) 


(6) 


With the composite Bragg theory and the wave spectrum, the backscatter from the ocean surface can be 
simulated for different incident angles and polarizations. 


2.2. Sea Surface Wave Model 
The variance spectrum for the Bragg scattering model is related to the directional wave number spectrum 
S({k,p) by: 


W(k, p)=S(k, p)+S(k, p+r) (7) 


Here we represent the wave spectrum with the saturation spectrum B (or the surface curvature spectrum), 
which is a function of the surface variance spectrum S: 


B(k) =K's(k) (8) 
where the rain-modulated sea surface wave spectrum B, ranges from short gravity waves to capillary waves 
is: 


Bs= Bgc_damped + Beap T Bring (9) 


and where Bgc damped is the gravity wave spectrum, damped by the precipitation, with a viscous dissipation 
rate for the rain-induced turbulence damping, Bcap is the spectrum in the gravity-capillary region, and Bring 
is the spectrum for the rain-induced ring waves. 


2.2.1. The Gravity and Gravity-Capillary Wave Spectra 
For the gravity and gravity-capillary region of the spectrum, the parameterization developed by Kudryavtsev 
et al. [2003] is: 


[Bot ks 0) + (B2(k, #) +4lpe(k, (2) 7)" ao 
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Beg(k, ọ)= 21/1 


where, k is the wave number, ¢ is the wind direction, n, « are two tuning parameters, |, is energy input due 
to generation of parasitic capillaries and is also a function of the effective growth rate [,(k,p), and where 
effectively, fy(k,ọp) is the difference between the wind growth rate £(k,p) and the viscous dissipation rate. In 
this model, the rain damping effect was added through the parameterization for the viscous dissipation 
rate. Therefore, the effective growth rate {,(k,0) may be described by: 


Ay,k? 
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(11) 
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where v, is the viscosity coefficient due to the rain damping, and we adopt the rain-induced damping 
model developed by Nystuen [1990]: 

vee (1 etd) TE + 9 (12) 
where ve is the rain-induced eddy viscosity. Tsimplis [1992] found v. to be about 3 X 10 ° m/s, when rain 
exists, no matter how high the rain rate is. Therefore, Contreras and Plant [2006] described ve with a function 
of the rain rate as: 
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Here / and v are differential length and velocity scales characteristic of the turbulence, C = 0.2095 is a con- 
stant coefficient, and /F is the energy flux into the eddy viscosity due to the rain impact: 
IF(D; RR) =N(D; RR) (= Pu D? w?(D)) W(D) (14) 


Here N is the rain drop size distribution (DSD), W is the terminal drop velocity, and D = I is the diameter of 
the rain droplet. 


For the gravity wave region, the two tuning parameters (n and «) are constants. With constant tuning 
parameters, Bgc damped Was Calculated by equation (10). However, for gravity-capillary waves, the tuning 
parameters are functions of wave number and then the Bap is computed, also using equation (10). Details 
of the two tuning parameters and intervals for the gravity and gravity-capillary wave numbers are the same 
as those of the KHCC spectrum [Kudryavtsev et al., 2003]. 

2.2.2. The Ring Wave Spectrum 

The ring wave spectrum as developed by Le Méhauté [1988] is also adopted in our model. Given a rain drop 
of radius R impinging upon water, the ring wave spectrum is: 


Bring (k, RR) = Sa [ N(D)W? (D) D? p S) a) 


0 


J1 Ce ag) 
8v, k? 


where q is the wave frequency, g is the gravitational acceleration, RR is the rain rate, D is the diameter of 
the raindrop, N(D) is the drop size distribution [Marshall and Palmer, 1948], W(D) is the terminal fall velocity 
of the drop, J2(kD/2) is the second order Bessel functions of the first kind, and the maximum limit of integra- 
tion t is set to 4. The ocean surface backscatter model including the simultaneous effects of wind and rain 
was developed using the sea surface spectrum in equation (9) and the composite Bragg theory. 

2.2.3. Results of the Sea Surface Wave Model 

Figure 2 shows the saturation spectra of wind waves at various wind speeds without the rain effect. Figure 3 com- 
pares the saturation spectra with the rain effect (solid lines), rain damping (dashed line), and the wind-wave-only 
spectrum (dotted line) for a wind speed of 10 m/s with different rain rates. The rain effect on the waves contains 
wind-generated waves damped by rain and the rain-generated ring waves. Regarding the magnitude of the 
spectral values, although rain-affected wave spectra are higher than the wave spectra only damped by rain (by 
adding the ring waves), the rain-affected wave spectra are still lower than the wind-induced wave spectra. Hence, 
the composite rain effect on the ocean surface waves is attenuation of these waves. 


2.3. Atmospheric Rain Effect on NRCS 

In the atmosphere, raindrops attenuate the radar transmitting signal and also induce volume scattering. 
These two mechanisms compete by decreasing or increasing the NRCS, respectively. An additive model, 
including the rain’s effect on the NRCS in the atmosphere is given by Nie and Long [2007]: 


Om =Osurf%atm + atm (16) 
-1 
10 
where Sm is the total NRCS received by 
radar, Gsurf is the backscatter from the 
>. a3 ms ocean surface, atm is the rain-induced 
g 10 E “13m/s atmospheric attenuation, and Gatm is 
= Sms the atmospheric volume backscatter. 
ba 
© 
£ 10°- 3 m/s nee 
3. Validation of the 
Composite Model 
10° 107 10° w To validate the accuracy of this wind- 
Wavenumber (rad/m) rain model, data sets of C-band VV- 


polarized NRCSs, wind speeds, and 
Figure 2. Model saturation spectra of short wind waves at various wind speeds i hürri 
without the rain effect: 3-23 m/s from lower to upper curves with the resolution rain rates over two hurricanes were 
of 5 m/s. investigated. 
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Figure 3. Saturation spectra for a wind speed 
trum with additive ring wave spectrum (solid). 


of 10 m/s with different rain rates: wind wave spectrum (dotted), damped wind wave spectrum (dashed), and damped wind wave spec- 


3.1. Wide-Swath SAR Data 

We have acquired two Envisat SAR images over Hurricane Gustav (03:56 UTC, 1 September 2008) and Hurri- 
cane Ike (04:23 UTC, 13 September 2008) from the European Space Agency (ESA). The SAR images are wide 
swath mode (WSM) with a medium resolution of 150 m and a swath width of 405 km at VV polarization. 
The ocean surface part of our model was based on Bragg resonance between the sea surface waves and 
the radiative waves; the sea surface waves were simulated by a wave spectrum, and therefore a particular 
scale of the ocean surface is needed to cover enough wave scales. Therefore, we calibrated the SAR image 
and then averaged the spatial resolution to 1 km (shown in Figure 4) with the boxcar averaging method. 
Moreover, with the averaging, the SAR speckle noise is almost removed. The calibration process is also a 
methodology to transfer the SAR image pixels to latitude and longitude spatial location information, using 
the transfer function given by ESA. 


3.2. SFMR-Measured Hurricane Wind and Rain 

At the two SAR imaging times, we also obtained wind speeds and rain rates, as measured by the stepped- 
frequency microwave radiometer (SFMR) on board the WP-3-D research aircraft, which is employed for 
NOAA's operational surface wind and rain measurements [Uh/horn and Black, 2003]. SFMR provides along- 
track mapping of wind speeds at relatively high spatial (1.5 km) and temporal (1 Hz) resolutions. These 
winds are well validated by measurements from dropwindsondes, each equipped with GPS. The RMS error 
for SFMR measurements is less than ~4 m/s at the surface and less than ~5 m/s at 10 m height [Uh/horn 
et al., 2007]. With six frequencies between 4.5 and 7.2 GHz used in the SFMR, the rain rate is also inferred 
with the wind speed. The rain rate error is less than ~5 mm/h, after algorithm correction and error analysis 
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Figure 4. Envisat SAR images for: (a) Hurricane Gustav (03:59 UTC, 1 September 
2008) and (b) Hurricane Ike (04:23 UTC, 13 September 2008). And the positions of 
SFMR adopted here were displayed in blue points as well as the SAR data 
matched up in cyan points. 


[iang et al, 2006]. Along the flight 
track (blue dots in Figure 4), longitude, 
latitude, sea surface wind speed, and 
rain rate are measured at the same 
time. 


3.3. Matchup SAR and SFMR for 
Moving Hurricanes 

SFMR data selected here were col- 
lected during a flight that lasted 1 h 
but two SAR image were captured 
almost instantaneously. During this 
1 h, the hurricane also moves. There- 
fore, our matchup process here was 
carried out to remove the effect of the 
relative movement between the air- 
craft carrying the SFMR and the hurri- 
cane. Moreover, we assume that the 
hurricane moves linearly along the 
Best Track (BT) during this 1 h. For the 
hurricane moving continuously along 
a trajectory, we map the hurricane 
center observed by SAR to the BT posi- 
tion at each SFMR measurement time. 
The center locations of BT and of SAR 
hurricane images detected by Li et al. 
[2013] were adopted, as shown in 
Table 1. To simplify the process, we 
assume that the hurricane moves line- 
arly between two BT locations without 
rotation. 


For each selected set of wind speeds 
and rain rates measured by SFMR, the 
mapping procedure consists of the fol- 
lowing steps: (1) interpolate the BT 


center locations to the time of the SFMR observation; (2) move the SAR image to colocate the SAR center 
and the interpolated BT center; (3) determine whether the SFMR data set is within the mapped SAR image, 
or not; (4) if yes, the data set from the SFMR and the NRCS for the same location in the mapped SAR image 
are considered to be matched (cyan dots shown in Figure 4). As shown in Figure 5, there are 296 matchups 
of wind speed, rain rate, and NRCS. The rain rates are between 10 and 35 mm/h, wind speeds are from 25 


to 45 m/s, and the NRCSs are between —10 and —5 dB. 


3.4. Validation Results 


To model the volume scattering and attenuation due to rain in the atmospheric column, three-dimensional 
rain rates are needed, which are difficult to measure. Therefore, we modeled the NRCSs for the matched 
data sets using only the ocean surface part. The model results with: (a) wind only (Figure 6a) and (b) with 


Table 1. Hurricane Center Information of SAR and Related Best Track Estimates 


SAR Best Track 
Name Date Time Latitude Longitude Time Latitude Longitude 
Gustav 1 Sep 2008 03:59 27.0°N 88.4°W 00:00 26.9°N 87.7°W 
06:00 27.9°N 89.0°W 
Ike 13 Sep 2008 04:23 28.0°N 94.6°W 00:00 28.3°N 94.0°W 
06:00 29.1°N 94.6°W 
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Z . f Rain 5 5 mm/hr both wind and rain (Figure 6b) are 
plotted against the matched data set. 
-6.5 For the wind-only model results, the 
Pa 30mm/hr comparison has a bias of 0.71 dB, a 
d] . $ standard deviation of 0.60 dB, and cor- 
Bigs 4 ” relation coefficient of 80.9%. When the 
S R 25 mmihr rain effect is taken into account, the 
N e 
2 8 o simulation results improved; the bias is 
z ; 0.20 dB with a standard deviation of 
<x 20mm/hr š : 
o -8.5) D 0.54 dB and correlation coefficient is 
eme 
° 81.9%. In the comparisons shown in 
-9} | : 
4Smamihr Figure 6, near the SAR NRCSs at about 
PP elles ta | —9.0 dB, the rain rate is about 
° 20 mm/h (hereafter defined as 
10 - - - 40mm/hr region-A), and our composite model 
25 30 35 40 45 - 
Wind speed (m/s) appears to simulate the attenuation 


of the NRCSs due to rain very well. 
When the NRCSs are about —7.7 dB, 
as observed by SAR, the effect of 
heavy rain is about 30 mm/h (hereafter defined as region-B), which also appears to be modeled appro- 
priately. However, when the SAR NRCSs is —7.0 dB with rain rate of 20 mm/h (defined as region-C), this 
case is not simulated very well, although the results presented in Figure 6b are still better than those of 
Figure 6a. In terms of winds, Figure 5 implies that the wind speeds are about 27.5 m/s for region-A and 
from 35.0 to 42.0 m/s for region-B, and resulting SAR NRCSs are reasonably well simulated. The disagree- 
ment between simulated and measured NRCSs occurs near wind speeds of about 32.0 m/s. This demon- 
strates that the bias of region-B is not caused by the model saturation due to high wind speeds. A 
possible reason may be the asymmetry and rotational motion of the hurricane. 


Figure 5. The collected matchup data sets: wind speed and rain rate measured 
by SFMR, as well as NRCS observed by SAR. 


4. Mechanisms for Rain Effects on SAR Hurricane Observations 


The modeling approach used here consists of modifying the parameterization for a wind wave spectra to 
include: (a) rain-induced turbulent damping and the enhancement of small gravity and gravity capillary 
waves due to the presence of ring waves, (b) the attenuation and volume scattering due to raindrops in the 
atmosphere, then (c) analyzing the backscatter caused by rain on both the ocean surface and in the atmos- 
phere. In total, there are four possible ways that rain can affect the NRCS: (1) rain droplets can attenuate the 
RADAR signal crossing the atmosphere; (2) rain droplets can cause volume scattering and thus increase the 
backscattering NRCS; (3) the radar waves can experience diffraction on rain products on the ocean surface 
including craters, crowns, and stalks; (4) rain changes the ocean surface waves as generated and driven by 
wind. In this paper, the former two ways are summarized as the atmosphere part, and the latter two ways, 
as the ocean surface part. 


4.1. Volume Scattering and Attenuation 
Based on the atmosphere model, the rain effect in the atmosphere is described by: 


Aco(dB)=10 : IG (surf %atm + atm) — 10 i 19 (S surf) (17) 


Therefore, C-band NRCS affected by rain in the atmosphere is simulated and shown in Figure 7 for the rain 
rates of 10, 30, and 50 mm/h, respectively, with the assumption that rain rates were not changing along the 
radiative incidence track. The VV-polarized hurricane measurements, which are mostly between —18 and 
—5 dB for winds above 5 m/s, are probably attenuated by the rain in the atmosphere; Figure 7 suggests 
that the effects of attenuation and volume backscattering in the atmosphere are less than 1 dB, even with a 
rain rate of 50 mm/h along the incidence track. As shown in Figures 8 and 9, when the NRCS is between 
—18 and —5 dB, the wind speed varies from 5 to 50 m/s. Moreover, the rain rate is not always the same 
along the entire radiative incidence track, as is often assumed in a hurricane; therefore, the real effect of 
rain in atmosphere should be less than the results simulated (with assuming a maximum rain rate of 
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Figure 6. Comparisons between simulated NRCS induced by the ocean surface 
part of the composite model and the ASAR observed NRCS: (a) simulation with 


only wind, and (b) simulation with wind and rain. 
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Figure 7. Rain effects on attenuation and volume backscattering in the atmos- 
phere as a function of the sea surface backscatter with three different rain rates. 


50 mm/h). Hence, we conclude that 
the atmospheric part can be neglected 
for C-band VV polarization in the hurri- 
cane observations. When the sea sur- 
face backscattering NRCS is less than 
—15 dB, representing wind speeds of 
about 40 m/s for the cross polarization 
(Figures10 and 11), the NRCS is 
enhanced because of the volume scat- 
tering, due to rain. 


4.2. Rain Products Diffraction 

For C-band W-polarized SAR observa- 
tions, the dependency of the NRCSs on 
the wind speed and the SAR geometries 
(e.g., incidence angle and azimuth angle 
with respect to wind direction) is gener- 
ally expressed by the GMF, such as the 
two-scale model [Wackerman et al., 2002] 
and CMODS for neutral winds (CMOD5.N) 
[Hersbach et al, 2007, Hersbach, 2010]. 
For cross-polarization images (VH, HV), an 
empirical function denoted C-2PO was 
developed by Zhang and Perrie [2012]. 
For the HH-polarization, no well- 
developed empirical function exists and 
the widely used approach is to develop a 
hybrid model function consisting of one 
of the GMFs and a polarization ratio [e.g., 
Zhang et al., 2011; Nunziata et al, 2014]. 
However, HH-polarization SAR images 
are not discussed in this study. 


We compare model results with the 
empirical functions, considering only 
the wind-induced KHCC wave spectrum 
for the VV-polarization SAR images and 
the cross-polarization SAR images. Fig- 
ures 8 (upwind) and 9 (crosswind) pres- 
ent the VV-polarized NRCSs simulated 
by this model, in comparison with 
results from CMODS5.N. One can see 
that the simulated NRCS results are 
close to those from CMODS.N for both 
the upwind and crosswind cases, illus- 
trating that this model can reliably 
simulate the NRCS, with respect to 
wind-generated KHCC wave spectra. 


However, for cross polarizations, we 
only simulate the modulation of NRCS 
by ocean wind without considering 
rain's effect. The wind-induced NRCS 
modulations in both upwind (Figure 
10) and crosswind (Figure 11) cases are 
much smaller than the results given by 
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Figure 8. NRCS as a function of the wind speed, comparing the composite model 
only, with KHCC wind wave spectrum (black circles) and the CMOD5N function. 
The incidence angle is 38°, and the relative wind direction is 0° (upwind) for 
C-band W-polarization. 


NRCS (dB) 
a 


O KHCC 
CMOD5N 
-25 1 4 4 4 L 
10 20 30 40 50 


Wind speed (m/s) 


Figure 9. NRCS as a function of the wind speed for the comparison between the 
composed model only with KHCC wind wave spectrum (black circles) and the 
CMODSN function. The incidence angle is 38°, and the relative wind direction is 
90° (crosswind) for C-band VV-polarization. 


4.3. Changes on the Wind-Induced Waves 


the C-2PO model. The reason is that 
the Bragg model is not suitable for 
wind-induced NRCS simulation for the 
cross-pol imagery, and certainly does 
not include the rain’s effect. Therefore, 
the rain's effect on ocean wave was 
not simulated for the cross-pol signal 
in this study. The underestimated 
NRCSs values may be caused by non- 
Bragg scattering, including specular 
reflection and diffraction of radio 
waves on sharp wedges of wave crest 
breaks [Kalmykov and Pustovoytenko, 
1976]. The non-Bragg scattering is not 
simulated in this model and is not 
related to polarizations. However, the 
VV-polarized NRCS values are much 
larger than those of cross-polarization 
images, with the same wind speed, 
which is evident, comparing Figures 8 
and 9 with Figures 10 and 11. Hence, 
although the non-Bragg scattering 
impacts on the cross polarization are 
significant, the impacts on the VV- 
polarization are slight. Just as for the 
crest breaks, we speculated that the 
diffraction of radio waves on the sharp 
wedges of craters, crowns, and stalks 
produced by rain colliding with the 
sea surface can be neglected for C- 
band VV polarized SAR hurricane 
observations, but should be important 
for the cross-polarization images. 


We have found that the previous three mechanisms for the rain effects on the NRCS can be neglected for 
C-band WV polarized SAR hurricane observations. Here we also consider the impacts of rain on wind- 
generated waves. Based on the validation in section 3.2, we simulated NRCSs with the rain damped wave 
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Figure 10. NRCS as a function of the wind speed for the comparison between the 
model only with KHCC wind wave spectrum (black circle) and the C-2PO function. 
The incidence angle is 38°, and the relative wind direction is 0° (upwind) for C- 
band cross polarization. 


spectrum, in comparison with the rain 
changed wave spectrum (that is, wave 
spectrum damped by rain and adding 
the ring waves generated by rain) for 
the C-band VV-polarization shown in 
Figure 12 for upwind wind cases, and 
Figure 13, for crosswind cases. There- 
fore, even with the rain-induced ring 
waves, the NRCSs are also attenuated 
by the precipitation, for both upwind 
(Figure12) and crosswind (Figure 13) 
directions. However, the attenuation 
of the NRCSs values due to the 
rain damping effect decreases, with 
increasing wind speeds. Conversely, 
the total rain effect on the sea surface 
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Figure 11. NRCS as a function of the wind speed for the comparison between the 
models only with KHCC wind wave spectrum (black circles) and the C-2PO func- 
tion. The incidence angle is 38°, and the relative wind direction is 90° (crosswind) 
for C-band cross polarization. 
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Figure 12. Simulated NRCSs with different wave spectra: only with wind wave 
spectrum (black circles), with the damped wind wave spectrum (blue points), with 
the damped wind wave and ring wave spectra (red circles), and the CMOD5.N 
function (solid line). The incidence angle is 38°, and the relative wind direction is 
0° (upwind) for C-band VV-polarization. From upper to lower limits on the plot, 
rain rates are from 5 to 60 mm/s with the increments of 5 mm/s. 
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Figure 13. Simulated NRCSs as Figure 12 but for crosswind (relative wind direc- 
tion is 90°). 


(wind wave damped by rain and the 
ring waves generated by rain) 
increases, with decreasing or during 
low wind speeds (<10 m/s), but 
decreases for increasing, moderate 
and high wind speeds. 


5. Conclusions 


In this paper, five possible mechanisms 
for the rain effects on the spaceborne 
C-band SAR hurricane observations are 
considered: (1) attenuation and (2) vol- 
ume backscattering for the microwave 
transfer in atmosphere, as well as (3) 
rain-induced damping to the wind 
waves and (4) rain-generated ring 
waves on the ocean surface, and (5) 
diffraction on the sharp edges of rain 
products. The first four mechanisms 
are simulated by a composite model, 
and the fifth mechanism was dis- 
cussed. Comparisons between the 
observed NRCSs and the atmosphere 
part of our model imply that two possi- 
ble effects in the atmosphere can be 
neglected for C-band VV-polarization 
SAR hurricane observations. However, 
the atmosphere part should be impor- 
tant at very low wind speed and the 
VV-polarized NRCS will be increased 
for the volume backscattering. Thus, 
we speculated that the diffraction of 
radio waves on craters, crowns, and 
stalks, induced by rain on the ocean 
surface can be neglected for C-band 
VV-polarization, and that the NRCSs 
are essentially affected only by wind 
and can be simulated well for C-band 
VWV-polarization data without including 
the diffraction of the radar by the 
sharp wedges of wave breaking but 
not well for the cross-polarization data 
(Figures 10 and 11). However, the 
mechanism of diffraction on the sharp 
wedges of wave breaking and rain 
products needs to be further modeled 
to validate this speculation. Our com- 
posite model was validated with the 
matchup of observations of wind 
speed and rain rate from the SFMR 
data as well as the NRCSs and inci- 
dence angles of C-band VV-polariza- 
tion SAR data over two hurricanes. 
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Therefore, the most important mechanism for the rain effect on the C-band VV polarized SAR hurricane 
observations is through the influence of waves on the ocean surface. Moreover, the non-Bragg parameter- 
ization will be added to our model to study whether or not it can explain the gaps in the C-2PO results, in 
future studies. 
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